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Finite-Length Multi-Input Multi-Output 
Channel Shortening Pre-filters 

Related Application: 

5 This application claims priority from Provisional application No. 60/158,713 

filed on October 8, 1999. This application is also related to a Provisional 
application No. 60/158,714, also filed on October 8, 1999. 

Background of the Invention: 

10 The combination of maximum likelihood sequence estimation (MLSE) with 

receiver diversity is an effective technique for achieving high performance over 
noisy, frequency-selective, fading channels impaired by co-channel interference. 
With the addition of transmitter diversity, the resulting multi-input multi-output 
':3 (MIMO) frequency-selective channel has a significantly higher capacity than its 
[ni5 single-input multi-output (SIMO) or single-input single-output (SISO) counterparts. 

The use of maximum likelihood multi-user detection techniques on these 
%3 frequency-selective MIMO channels significantly outperforms single-user detection 
1" techniques that treat signals from other users as colored noise. However, MLSE 

complexity increases exponentially with the number of inputs (or transmit 
ILSO antennas) and with the memory of the MIMO channel, making its implementation 
J:^ over sever inter-symbol interference (ISI) channels very costly, 
u The MIMO channel can be modeled as a collection of FIR filters (i.e., an 

FIR filter between each input point (e.g., transmitting antenna) and each receiving 
point (e.g. receiving antenna), and the "memory of the channel" corresponds to the 
25 number of taps in the FIR filters. 

The Discrete Matrix Multitone (DMMT) was shown to be a practical 
transceiver structure that asymptotically achieves the MIMO channel capacity 
when combined with powerful codes. It uses the Discrete Fourier Transform 
(DFT) to partition the frequency responses of the underlying frequency-selective 
30 channels of the MIMO systems into a large number of parallel, independent, and 
(approximately) memoryless frequency subchannels. To eliminate inter-block and 
intra-block interference, a cyclic prefix whose length is equal to the MIMO channel 
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memory is inserted in every block. On severe-ISI MIMO channels, the cyclic prefix 
overhead reduces the achievable DMMT throughput significantly, unless a large 
FFT size is used which, in turn, increases the computational complexity, 
processing delay, and memory requirements in the receiver. 

In short, the computation complexity increases exponentially with the 
number of taps in the FIR filters that may be used to model the channel. 

N. Al-Dhahir and J. M. Cioffi, in "Efficiently-Computed Reduced-Parameter 
Input-Aided MMSE Equalizers for ML Detection: A Unified Approach," IEEE Trans. 
Information Theory, pp. 903-915, May 1996, disclose use of a time-domain pre- 
filter in the receiver to shorten the effective channel memory and hence reduce the 
cyclic prefix overhead and/or the number of MLSE states. The disclosed 
approach, however, is for SISO systems, and not for MIMO systems. 

Summary 

An advance in the art is realized with a MIMO pre-filter that is constructed 
from FIR filters with coefficients that are computed based on environment 
parameters that are designer-chosen. Given a transmission channel that is 
modeled as a set of FIR filters with memory v , a matrix W is computed for a pre- 
filter that results in an effective transmission channel B with memory A//,, where 
A/^ < V , where B is optimized so that 8^,^^ = argmin^ traceiR^^) subject to selected 

constraints; Ree being the error autocorrelation function. The coefficients of W, 
which are sensitive to a variety of designer constraints, are computed by a 
processor within pre-filter at the front end of a receiver and loaded into an array of 
FIR filters that form the pre-filter. 

Brief Description of the Drawing 

FIG. 1 shows the major elements of a receiver in accord with the principles 
disclosed herein; 

FIG. 2 presents the structure of pre-filter 30; and 

FIG. 3 is a flowchart describing the method carried out by processor 220 
within pre-filter 30. 
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Detailed Description 

FIG. 1 depicts the general case of an arrangement with f7q transmitting 
antennas 11-1, 11-2, ... Il-Hq, that output signals (e.g., space-time encoded 
5 signals) to a transmission channel, and Dq receiving antennas 21-1, 21-2, ... 21 -/Iq. 
Each transmitting antenna p outputs a complex-valued signal Xp, the signals of the 
Hq antennas pass through a noisy transmission channel, and the no receiving 
antennas capture the signals that passed through the transmission channel. The 
received signals are oversampled by a factor of / in element 20 and applied to pre- 
10 filter 30. Thus, the sampling clock at the output of element 20 is of period Ts-T/I, 
where 7 is the inter-symbol period at the transmitting antennas. The transmission 
channel's characterization is also referenced to Tg. In the illustrative embodiment 
^3 disclosed herein, therefore, pre-filter 20 develops n\ output signals that are applied 
£n to a conventional multi-input receiver 40, and the received signal can be 
;'=Ji5 expressed by 

5 yl^^ = ii;e'xi!^.+ni>'. (1) 

Q where y]J^ is the signal at time k at the receiving antenna, Z?^'^' is the m*^ 
=rj coefficient (tap) in the channel impulse response between the transmitting 

[y antenna and the receiving antenna, and n^^^ is the noise vector at the 

11 

[20 receiving antenna. The memory of this path (i.e., the largest value of m for which 
hli;^^ is not zero) is denoted by v^''^^ . It not unreasonable to assume, however, that 
the memory of the transmission channel is the same for all ij links {n^xn^ such 

links), in which case v^''^^ = v . Alternatively, the v^''^^ limit in equation (1) can be 
set to that v which corresponds to maximum length of all of the x channel 

25 input responses, i.e., v = max^y v^'*^* . All of these variables in equation (1) are 

actually /x1 column vectors, corresponding to the / time samples per symbol in 
the oversampled FIG. 1 arrangement. By grouping the received samples from all 
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Ho antennas at symbol time k into an nJx^ column vector , one can relate 
to the corresponding n, x 1 (column) vector of input samples as follows 

V 

y* = S"mX,-„+n,. (2) 

where Hm is the MIMO channel coefficients matrix of size nj x nj , x^_^ is a size 
n, x1 input (column) vector, and n^^ is a size njxl vector. 

Over a block of Nf symbol periods, equation (2) can be expressed in matrix 
notation as follows: 
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or, more compactly, 

Vk+Nf-r.k ~ 



•k+Nf-r.k-v '^^k+Nf-tk ' 



(4) 



The subscripts in equation (4) indicate a range. For example k-\-Nf~1:k 
indicates the range k^N^-^ and /c, inclusive. 

It is useful to define the following correlation matrices: 

^xy = ^{^k+Nf-tk-vY k^Nf-\M \ " ^xx^ 

^yy = ^[yk+Nf-tkYk+Nf-A-.k] ~ '^^xx'^ + ^nn ' 



(5) 
(6) 
(7) 
(8) 



^xx — ^[^k+Nf-U-v^k+Nf-tk-v] 
^nn = ^[^k+Nf-tk^k+Nf-tk] • 

It is assumed that these correlation matrices do not change significantly in time or, 
at least, do not change significantly over a time interval that corresponds to a 
TDMA burst (assumed to be much shorter than the channel coherence time), 
which is much longer than the length of the pre-filter, in symbol periods denoted by 
A/f. Accordingly, a re-computation of the above matrices, and the other parameters 
disclosed herein, leading to the computation of pre-filter coefficients, need not take 
place more often than once every TDMA burst. 
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Once H, Rxx and Rnn and known, Rxy and Ryy are computed by R^^H* and 
HR,X+R,„, respectively. 

Given the MIMO channel matrix H with v + 1 members (Ho,Hi,"-HJ, the 
objective is to create a MIMO pre-fiiter W (element 30 in FIG. 1) with M matrix 
5 taps, i.e., matrix W = [Wq W, W/y^.^]'^ , that equalizes H so as to create an 

overall transmission channel for receiver 40 that corresponds to a matrix B with 
memory , where N^«y . 

The matrix B can be expressed as B = [Bq B^ • • B;^^ ]^ where each Bj 
is of size x n^. . 

10 The MIMO channel-shortening pre-filter W (element 30) is conditioned, or 

^4 adjusted, to minimize the equalization Mean Squared Error (MSE), defined by 

In MSE = ^race(Reg) , where Ree is the autocorrelation matrix of the error vector Ek 

in 

iQ that is given by 

%j E^ = B x^^yy^_^.^_^ - W )ff^^[^^_^.i^ , (9) 

s 15 where the augmented MIMO matrix, B* , is 

B^-[On.x..A b; b; ... b;;;^ o,^,„^,]^[o,^,,^, b* o„^,,^,], (io) 

f U A is the decision delay that lies in the range 0 < (A/^ + v - A/^ - 1) , and 
Q s = Nf+v-N^-A-^. The /?, x /?, error autocorrelation function Ree can be 
expressed by the following: 

= B*R^B 
= B'RB, 

where R is a sub-matrix of R^ determined by A . 

Using the orthogonality principle, which states that E[E;,y^^„^_^^] = 0 it can 

be shown that the optimum channel-shortening pre-filter and target impulse 
response filters (W and B, respectively) are related by 
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= B;,,R,,H(HR.X+Rnnr (12) 

The last line shows explicitly that the MIMO channel-shortening pre-filter consists 
of a noise whitener R~l , a MIMO matched filter H*, and a bank of FIR channel- 
shortening pre-filter elements. 

5 It remains to optimize B such that the MSE is minimized, which may be 

obtained by computing the parameters of B that, responsive to specified 
conditions, minimizes the trace (or determinant) of Ree- The following discloses 
two approaches to such optimization. 

Under one optimization approach the coefficients of B are constrained so 
Cao that some coefficient of B is equal to the identity matrix, I. A solution subject to 
f = this Identity Tap Constraint (ITC) can be expressed by 
;;i BJJ = argming frace(RgJ subject to 6*0) = l„ , (13) 

J: where <P* = [0,^,„^„ On,x„,(A,,-m)] and 0 < m < A/, . It can be shown that the 

l"^ optinnum MIMO target impulse response and the corresponding error 

^^15 autocorrelation matrix are given by 

rU B[J = R-'<l>(0*R-'a))-' and (14) 

" R':^^^={O'R-'0r\ (15) 



As indicated above, R is affected by the delay parameter A. Unless dictated by 

the designer, the delay parameter A , which can range between 0 and 
20 (A/^ + V - A/^ - 1) , is chosen to minimize the trace of R'ee^n • Similarly, the index 

parameter m, which ranges between 0 and N^^, and which that affects <P , is 

chosen to minimize the trace of Ree^^nin ■ 

Under a second optimization approach the imposed constraint is B*B = 1^^ . 

A solution subject to this Ortho-Normality Constraint (ONC) can be expressed by 
25 B^pf^ = argming trace{R^^) subject to B*B = 1^^ , (16) 

Defining the eigen-decomposition 
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R^UZU* =Ud/ag(<7o.cr,-a„,^^,,.,)U-. (17) 
where ctq > •• > Cn.tw.+iH . then the optimum MIMO target response and the 
resulting error autocorrelation matrix are given by 

B^;;^=U[e,,^ ... e,,^^,,,]. (18) 
5 where e, is unit vector with a 1 at position /, and O's elsewhere, and 

RSn =c'/agKA/,^■^^T,^(A/,.1)-1^ (19) 
Illustratively, if a7/=3 and A/b=3, B^^^ = [J[eg,e^o,e^^], meaning that is a three 

column matrix comprising the 9**^ through the 1 1*^ columns of matrix U. Stated in 
words, the optimum MIMO target impulse response matrix is given by the n\ 
10 eigenvectors of R that correspond to its n\ smallest eigenvalues. The delay 
parameter A (0 < A < A/^ + v - A/^ - 1) that affects R is optimized to minimize the 

trace (or determinant) of Ree% • 

With the above analysis in mind, a design of a prefilter 30 can proceed for 
any given set of system parameters, which includes: 
Cas • MIMO channel memory between the input points and the output point of the 
□ actual transmission channel, v , 

The number of pre-filter taps chosen, A/f, 
The shortened MIMO memory, A/b, 
The number of inputs to the transmission channel, Hj, 
The number of output derived from the transmission channel, /?o, 
The autocorrelation matrix of the inputs, Rxx, 
The autocorrelation matrix of the noise, Rnn, 
The oversampling used, /, and 
The decision delay, A . 
25 The structure of pre-filter 30 is shown in FIG. 2. In the illustrated 

embodiment, it comprises two main components: processor 220 and filter section 
210. 

Filter section 210 in the FIG. 2 illustrative example comprises a collection of 
FIR filters that connect the rio input array of signals from sampling circuit 20 to an 
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Hi output array of points. That is, there are ) x / FIR filters Py, , , that couple input 
point y to output point /. 

Processor 220 is responsive to the no signals received by antennas 21 and 
sampled by circuit 20, and it computes the coefficients of W, as disclosed above. 
Wo is a matrix that defines the coefficients in the 0*^ tap of the j xi FIR filters, Wi 
is a matrix that defines the coefficients in the 1®* tap of the j xi FIR filters, etc. 

The method of developing the parameters of pre-filter 30, carried out in 
processor 220, is shown in FIG. 3. Block 100 develops an estimate of the MIMO 
channel between the input points and the output point of the actual transmission 
channel. This is accomplished in a conventional manner through the use of 
training sequences. The estimate of the MIMO channel can be chosen to be 
limited to a given memory length, v , or can be allowed to include as much 
memory as necessary to reach a selected estimate error level. That, in turn, 
depends on the environment and is basically equal to the delay spread divided by 
Ts. 

Following step 100, step 110 determines the matrices, Rnn, Rxx, Rxy, and 
Ryy. The matrix Rnn is computed by first computing n = y - Hx and then computing 

the expected value E[n*n] -- see equation (8) above. The matrix Rxx is computed 
from the known training sequences see equation (7) above (or is pre- 
computed and installed in processor 220). In may be noted that for uncorrelated 
inputs, Rxx=l. The matrices Rxy and Ryy are computed from the known training 
sequences and the received signal or directly from H and Rnn- see equations (5) 
and (6) above. 

Following step 110, step 120 computes R-^ = R^^, -R^^R^R^ , and the sub- 
matrix R . From equation (10) is can be seen that R is obtained by dropping the 
first A rows and the last n,s rows of R^ . 

In accordance with the ITC approach, selecting some value of 0 < m < A/^ 
allows completion of the design process. Accordingly, following step 120, step 
130 chooses a value for m, develops O* = [0„^^„^^ 1^^ ^n,>cn,(A/^-m)] carries out 
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the computation of equation (13). Step 140 finally develops the coefficients of 
matrix W in accordance with equation (12), and installs the developed coefficients 
within filter 210. 

In accordance with the ONC approach, step 130 computes the matrix U in a 
5 conventional manner, identifies the unit vectors e/, and thus obtains the matrix B. 
Step As with the ITC approach, step 140 develops the coefficients of matrix W in 
accordance with equation (12), and installs the developed coefficients within filter 
210. 

It should be understood that a number of aspects of the above disclosure, 
10 for example, those related to the ITC constraint and to the ONC constraint, are 
merely illustrative, and that persons skilled in the art may make various 
modifications that, nevertheless, are within the spirit and scope of this invention. 
^ For example, the pre-filter described above generates a multi-output signal, with 
in the number of outputs being /?/, that being also the number of transmitting 
fds antennas 11. This, however, is not a limitation of the principles disclosed herein. 
'^^ The number of pre-filter outputs can, for example, be larger than n,, for example as 
Q high as n, (A/^ +1) . The performance of the receiver will be better with more filter 

Q outputs, but more outputs require more FIR filters, more FIR filter coefficients, and 
Jr? correspondingly, a greater processing power requirement placed on processor 
fi^O 220. 
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